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ABSTRACT 

We present the first results from a high samphng rate, multi-month rever- 
beration mapping campaign undertaken primarily at MDM Observatory with 
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supporting observations from telescopes around the world. The primary goal 
of this campaign was to obtain either new or improved H/5 reverberation lag 
measurements for several relatively low luminosity AGNs. We feature results 
for NGC 4051 here because, until now, this object has been a significant out- 
lier from AGN scaling relationships, e.g., it was previously a ~2-3(T outlier on 
the relationship between the broad-line region (BLR) radius and the optical 
continuum luminosity — the i?BLR^-^ relationship. Our new measurements of 
the lag time between variations in the continuum and H/3 emission line made 
from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of 
i?BLR = l.S7toit light days and black hole mass of Mbh = (1.73lg:|) x IO^Mq. 
This radius is consistent with that expected from the Rblr^L relationship, based 
on the pres ent luminosi t y of N GC 4051 and the most current calibration of the 



relation by iBentz et al.l (l2009al ). We also present a preliminary look at velocity- 
resolved H/3 light curves and time delay measurements, although we are unable 
to reconstruct an unambiguous velocity-resolved reverberation signal. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 



Recent theoretical and observational studies have provided strong evidence suggesting 
a connection between supermassive black hole (SMBH) growth and galaxy evolution (e.g ., 



Bennert et al. 2008: Somerville et al. 2008: Shankar et al. 



20091 : iHopkins fc HernquistI 120091 ). 



To better understand this connection, we need more direct measurements of SMBH masses 
across cosmological distances. Unfortunately, measuring SMBH masses directly with dy- 
namical methods requires high angular resolution, so use of these methods is limited to 
relatively nearby sources. This resolution problem can be obviated by studying SMBHs in 
type 1 active galactic nuclei (AGNs) . In t his case, the technique of reverberation mapping 
(IBlandford fc McKedll982l : |Petersonlll993l ) may be applied to measure the SMBH mass, as 
has been done for ra ore than 3 dozen type 1 AGNs to date (e.g., see recent compilation by 
Peterson et al1l2004l ). 



Reverberation mapping relies on time resolution rather than angular resolution, since it 
takes advantage of the presence of a time delay, r, between continuum and emission line flux 
variations observed through spectroscopic monitoring. This time delay corresponds to the 
light travel time across the broad-line region (BLR), and thus measurements of r provide 
an estimate of the size of the region, -Rblr = cr. Because the BLR gas is well within the 
sphere of influence of the black hole and studies have provided evidence for virialized motions 
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within this region (e.g.. iPeterson et al.ll2004l . and references therein), -Rblr can be related to 
the mass of the SMBH through the velocity dispersion of the BLR gas. 

Although stellar and gas dynamical modeling and reverberation mapping are invaluable 
for measuring SMBH masses directly, these methods are observationally taxing, due to reso- 
lution requirements for dynamical methods and time requirements for reverberation mapping 
campaigns. It is currently impossible to measure SMBH masses directly for large statistical 
samples of galaxies. However, these direct mass measurements have led to the discovery 
of scaling relationships that relate SMBH mass to other galaxy or AGN observables that 
can be used to investigate the connection between SMBH mass and galaxy evolution. In 
particular, some relations show connections between properties of the SMBH (i.e., its mass) 
and global properties of the host galaxy. Examples include the correlation between SMBH 
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tions connect various AGN propert ies. One such relation is the correlation between black 



hole mass and optical luminosity (iKaspi et al.l l2000l : iPeterson et al.l |2004| ). which relates 



directly to the accretion rates of AGNs. There is also a correlat i on between BLR radius 



and A GN luminosity, i.e., the Rblr^L relation (IKaspi et al.ll2000l . l2005l : iBentz et al.l 12006 



2009al ). which has proven to b e very powerful for mak i ng indirect SMBH mass estimations 
from single-epoch spectra (e.g . , Vestergaard I2OO2I . [2004 Corbett et al. 2003 ; KoUmeier et al. 



20061 : IVestergaard et all I2OO8I : IShen et al 



2008 



iFine et al.l 



2008h . 



These indirect mass 



estimates can then be related to other properties of the host galaxy through direct measure- 
ments or separate scaling relations. 

Although scaling relations have become widely used for statistical studies, it is important 
to understand that the indirect mass estimates determined by these relations are only as 
reliable as the direct mass measurements used to calibrate them. Therefore, establishing a 
secure calibration across a wide dynamic range in parameter space and better understanding 
any intrinsic scatter in these relations is essential. To accomplish this, we must continue to 
make new direct measurements as well as to check previous results that are, for one reason 
or another, suspect. 

NGC 4051, an SAB(rs)bc galaxy with a narrow-line Seyfert 1 (NLSl) nucleus at redshift 
z = 0.00234, is a case in p oint. Measurements of the BLR radius and optical luminosity 
(IPeterson et al.ll2000l . 120041 ) place it abo ve the -Rrt,r-L re lation, i.e., the BLR radius is too 
large for its luminosity (cf. Figure 2 of iKaspi et al.ll2005l ). It also appears to be accreting 
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mass at a lower Eddington rate than other NLSls (cf. Figure 16 of iPeterson et aLll2004j ). 
Th ese two anomal i es tog e ther s uggest that perhaps the BLR radius has been overestimated 
by iPeterson et al.l (l200d . 120041): indeed an iiidepend ent reverberation measurement of the 



BLR radius i n NGC 4051 by 



Shemmer et al.l (l2003l . hereafter 803) is about half the value 



measured by IPeterson et al. J2000l . hereafter POO). Furthermore neither the POO nor 803 



data sets are particularly well sampled on short time scale s, so ne i ther s et is suitable for 
detection of smaller time lags (e.g., <2-3 days). In addition. iRusselll (120031 ) reports a TuUy- 
Fisher distance to NGC 4051 that is ~50% larger than that inferred from its redshift (i.e., 
15.2 Mpc versus 10.0 Mpc, respectively). This suggests that the luminosity derived in past 
studies from the redshift could be an underestimate and might also be a contributing factor 
to the placement of NGC 4051 above the Rblr^L relation. 

In this work, we present an analysis of new, optical spectroscopic and photometric 
observations of NGC 4051, which represent the first results from a densely sampled rever- 
beration mapping campaign that began in early 2007. The campaign spanned more than 4 
months, during which time we consistently obtained multiple photometric observations per 
night and spectroscopic observations nearly every night from a combination of five different 
observatories around the globe. The immediate goal of this campaign is to remeasure the 
H/3 reverberation lag measurements for several objects on the low-luminosity end of i?BLR^-^ 
scaling relationship with poorly determined reverberation lags and, consequently, poorly de- 
termined black hole masses. We will also add to the overall sample of reverberation mapped 
AGNs by measuring lags for new objects. Another goal for this extensive collection of homo- 
geneous data is to reconstruct the observed response of the H/3 emission line to an outburst 
from the variable continuum source by modeling the response as a function of bot h line-of- 
sight velocity and light travel time, i.e., a "velocity-delay map" (for a tutorial, see IPeterson 
200ll : iHorne et al.ll2004l ). Creation of a velocity-delay map will provide novel insight into 
the structure and kinematics of the BLR. Though we have not yet attempted to reconstruct 
a full velocity-delay map, we present preliminary velocity-resolved lag measurements for 
NGC 4051. Complete results for NGC 4051 and other campaign targets will be presented in 
future work. 



2. Observations and Data Analysis 



Most data acquisit ion and analysis practices employed he re follow closely those described 
by lDenney et al.l (120061 ) and laid out by IPeterson et al.l (120041 ). The reader is referred to these 
works for additional details and discussions. Throughout this work, we assume the following 
cosmological parameters: Qm = 0.3, = 0.70, and Hq = 70 km sec~^ Mpc~^. 
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2.1. Spectroscopy 

Spectra of the nuclear region of NGC 4051 were obtained from both the 1.3-meter tele- 
scope at MDM Observatory and the 2.6-meter Shajn telescope of the Crimean Astrophysical 
Observatory (CrAO). The MDM observations utilized the Boiler and Chivens CCD spectro- 
graph, where 86 observations were taken over the course of 89 nights between JD2454184 
and JD2454269, targeting the H/?A4861 and [Olll] AA4959, 5007 emission line region of the 
optical spectrum. The position angle was set to 0°, with a slit width of 5'.'0 projected on 
the sky, resulting in a spectral resolution of 7.6 A across this emission- line region. Figure 
[U shows the mean and rms spectra of NGC 4051 based on the MDM observations. We 
acquired 22 CrAO spectra over 34 nights between JD2454266 and JD2454300 with the Na- 
smith spectrograph and SPEC-10 1340 x 100 pixel CCD. For these observations a 3'.'0 slit 
was utilized, with a 90° position angle. Spectral wavelength coverage for this data set was 
from ~3800-6000 A, with a dispersion of 1.8 A/pix and a spectral resolution of 7.5 A. Note 
that (1) the dispersion varies with wavelength: the value 1.8 A/pix is given for 5100A, and 
(2) the real wavelength coverage is slightly greater than given but the red and blue edges of 
the CCD frame are unusable (too low S/N ratio) because of vignetting. 

A relative flux calibration of each set of spectra was performed based on the constancy of 
the narrow [O III] A5007 line flux. Because this line emission originates in the extended, low- 
density narrow line region, it can be assumed constant over the timescale of this campaign and 
therefore serves as the basis for a reliable relative flux calibration. However, the data quality 
is not identical from night to night due to, e.g., seeing, weather conditions, atmospheric 
transparency, etc. This affects not only the integrated line flux in each observation, but also 
properties of the spe ctrum such as S/N and resolution. Consequently, we employ the spectral 



scaling algorithm of Ivan Groningen &: WandersI (Il992l ) for the [O ill] A5007 flux calibration. 



This algorithm determines the best scaling through minimization of residuals rather than 
simply calculating a simple multiplicative scale factor to scale the spectral fluxes. Following 
this method, we created a reference spectrum by averaging all spectra from a given data 
set. We then formed a difference spectrum by subtracting the reference spectrum from each 
individual spectrum. The algorithm uses a least squares method to minimize the residuals of 
the [O III] A5007 line flux in each difference spectrum by making small zero-point wavelength 
calibration adjustments, correcting for resolution differences, and applying a multiplicative 
scale factor to the [O III] A5007 line flux of the individual spectrum. Because this method is 
based on minimizing residuals between each individual spectrum and the reference spectrum, 
there is a small residual dispersion in the line fluxes after calibration. This dispersion is 
related to the data quality and the ability of the scaling algorithm to mitigate night to night 
differences b etween individual spectra, related largely to seeing and S/N. Tests of the original 



algorithm by Ivan Groningen fc WandersI (119921 ) estimated errors in the scaled fluxes of better 
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than 5%, however, past studies employing somewhat improved versio ns of this same scahn g 



algorithm typically achieved dispersions across a data set of ~2% (e.g. JPeterson et al.lll995l ). 
We measure a dispersion of ~1.5%, demonstrating the high level of homogeneity that we 
have been able to achieve in the current data set. 



2.2. Photometry 



In addition to spectral observations, we obtained supplemental \^-band photometry 
from the 2.0-m Multicolor Active Galactic NUclei Monitoring (MAGNUM) telescope at the 
Haleakala Observatories in Hawaii, the 70-cm telescope of the CrAO, and the 0.4-m telescope 
of the University of Nebraska. 

The MAGNUM o bservations were imaged w ith the multi color imaging pho t omete r 
(MIP) as described bv iKobavashi et al.l Jl998al lb[). lYoshiil J^OoJ, and IVoshii etaP J2003h . 
Photometric fluxes measured from 23 observations between JD2454182 and JD2454311 
within an aperture of 8'.'3. Photometric reduction of NGC 4051 was similar to that de- 



scribed for other sources by lMinezaki et al.l (120041 ) and ISuganuma et al.l (120061 ) . except the 
host-galaxy contribution to the flux within the aperture was not subtracted and the filter 
color term was not corrected because these photometric data were later scaled to the MDM 
continuum light curve (as described below). Also, minor corrections (of order 0.01 mag or 
less) due to the seeing dependence of the host-galaxy flux were ignored. 

The 76 CrAO photometric observations were collected between J245D4180 and JD2454299 
with the AP7p CCD mounted at the prime focus of the 70-cm telescope (/ = 282 cm). In 
this setup, the 512 x 512 pixels of the CCD field covers a 15'. x 15'. field of view. Photometric 
fluxes were measured within an aperture of 15'.'0. For further d etails of the CrAO V ^-band 
observations and reduction, see the similar analysis described by ISergeev et al.l (120051 ). 



The University of Nebraska observations were conducted by taking and separately mea- 
suring a large number of one-minute images (~20) each of 28 nights between JD2454195 
and JD2454290. D etails of the observing and reduction procedure are as described by 



Klimek et al.l (1200411. Comparison star m agnitudes were calibrated following lDoroshenko et al. 
(2005a 3) and Chonis fc Gaskell (2008). To minimize the effects of variations in the image 



quality, fluxes were measured through an aperture of radius 8 arcseconds. The errors given 
for each night are the errors in the means. 
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2.3. Light Curves 

Light curves of tlie H/5 flux were made based on integrated fluxes measured in the MDM 
and CrAO spectra between 4815-4920 A and over a linearly interpolated continuum deflned 
between the average flux density in each of the following regions blueward and redward of H/5, 
respectively: 4770-4780 A and 5090-5130 A. The CrAO H/? light curve was then scaled to 
the MDM light curve with a multiplicative constant based on the average flux ratio between 
the four pairs of closely spaced points in the MDM and CrAO H/? light curves separated by 
no more than 0.5 day. This scaling is necessary to account for differences in the amount of 
[O III] A5007 emission line flux that enters the slit in the different data sets due to seeing and 
aperture affects. The lower panel of Figure [2] shows the H/3 light curve derived from both 
data sets after scaling the CrAO fluxes to those measured from MDM spectra. 

A continuum light curve was created with observations from each l^-band photometric 
data set and the average continuum flux density measured over 5090-5130 A (i.e., rest 
frame ~5100 A) in each spectrum of the spectroscopic data sets. First, the multiplicative 
scale factor determined for above to scale the CrAO H/3 fluxes to the MDM light curve 
was also applied to the CrAO continuum fluxes, since the calibration of these fluxes is also 
susceptible to the same seeing and aperture affects as the H/? flux calibration. Next, this 
light curve as well as the individual photometric light curves (see Fig. ^ upper panel) 
were scaled, one by one, to the same flux scale as the MDM light curve by making an 
additive, relative flux adjustment to each. This additive correction is necessary for both 
spectroscopic and photometric data because of differences in host galaxy starlight that enters 
the different aperture sizes of the various data sets. For the photometric observations, there 
is an additional component (also additive) due to the larger width of the fllter bandpass. 
Each light curve was merged with the parent light curve to which it was scaled before the 
next light curve was scaled, thus building up a larger, more well-sampled light curve in the 
following order: MDM, MAGNUM, CrAO photometry. University of Nebraska, and CrAO 
spectroscopy. This was done so that the smaller, and in some cases shorter, light curves 
could be scaled to a longer and more densely sampled parent light curve. The scale factor 
applied to each secondary light curve to scale it to its parent light curve was calculated based 
on the difference between a linear least squares flt to this light curve and to the parent light 
curve before it, starting with the MDM light curve as the initial parent light curve. The 
flts to each light curve, both parent and secondary were limited to using only observations 
within the same overall temporal range, so that, when necessary, the beginning and/or ends 
of the light curves were truncated during the fltting process. 

Assignment of uncertainties to the photometric fluxes is described above in the text or 
in references describing the photometric data sets. However, we calculated the uncertainties 
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in the MDM and CrAO spectroscopic fluxes after creating tlie liglit curves for these data 
sets but before intercahbration. Typically, we determine uncertainties in our light curve 
flux measurements by applying a mean fractional error to all points. This fractional error is 
determined by comparing the average flux difference between closely spaced pairs of obser- 
vations, assuming that flux differences across these short times scales ar e due to noise rather 



than genuine variability. Because real variability has been confirmed by lKlimek et al.l (120041 ) 
to occur in NGC 4051 on time scales shorter than 2 days, and our sampling rate is ~1 day, 
on average, we could not use this method to determine the relative errors on our spectro- 
scopic flux measurements for this object. Instead, we took advantage of the observations of 
other higher-luminosity AGNs that we monitored as part of this larger campaign (i.e., same 
telescopes, instrumental setup, and observing conditions; results in preparation). Unlike 
NGC 4051, these objects neither exhibit variability on such short time scales nor have such 
short measured lags. Therefore, the uncertainties assigned to the NGC 4051 spectroscopic 
observations seen in Figure [2] are an average of the uncertainties in the flux measurements 
calculated as described above from closely spaced observations (separations of <2.0 days) of 
these other objects (e.g., typical fractional errors in the range ~ 0.12 — 0.21 and ~ 0.14 — 0.27 
for the continuum and line fluxes, respectively). 

The merged continuum and H/5 light curves shown in Figure [3] are used for the subse- 
quent time series analysis. These differ from simply combining the individual light curves, 
shown in Figure M, in the following ways: 

• First, we applied an absolute flux calibration to both light curves by applying a single 
multiplicative scale factor determined from the ratio of the [O ill] A5007 emission line 
flux determined by POO to that measured in the reference spectrum used above for the 
relative flux calibration. Unlike the emission line flux in our reference spectrum, the 
POO flux measurement of F([0 III] A5007) = (3.91 ± 0.12) x 10"^^ erg s"^ cm'^ was 
taken from observations obtained under photometric conditions and, consequently, 
~8% larger than our measured value. Additionally, the POO measurement was made 
employing observing strategies and measurement practices similar to what we present 
here, thus validating this direct comparison. This additional flux calibration does not 
affect the reverberation results but is necessary for accurately measuring the 5100A 
continuum luminosity. 

• Second, we subtracted the host g alaxy starligh t contri bution to the continuum flux. 



determined using the methods of iBentz et al.l (l2009al ) to be Fgai{5100A)= (9.18 ± 
0.85) X 10"^^ erg s"^ cm-^ A-\ 

Third, we binned closely spaced observations as a weighted average and applied this 
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to continuum flux measurements separated by <0.25 days and H/5 flux measurements 
separated by <0.5 days. 



Fluxes for individual observations (i.e., before time binning) from all sources are listed 
in Table [TJ Values listed represent the flux of each observation after completing all flux 
calibrations described above (i.e., relative calibration to intercalibrate all data sets onto 
the MDM flux scale, followed by absolute calibration based on the POO [O ill] A5007 line 
flux and removal of host starlight contamination). Column 1 gives the Julian Date of each 
observation. The 5100 A continuum or V-band flux and integrated H/5 flux are given in 
columns (2) and (3), respectively, and column (4) lists the source of each measurement. 
Photometric and spectroscopic observations from CrAO can be differentiated by noting that 
no B.(3 flux values are present for photometric observations. 

Table [2] displays statistical parameters describing the final light curves shown in Figure 
El Column (1) gives the spectral feature represented by each light curve, and the number of 
data points in each light curve is shown in column (2). Columns (3) and (4) are mean and 
median sampling intervals, respectively, between data points. The mean flux with standard 
deviation is given in column (5), while column (6) shows the mean fractional error in these 
fluxes. Column (7) gives the excess variance, calculated as 



{/) " 

where is the variance of the observ ed fluxes, 5'^ is their mean square uncertainty, an d 
(/) is the mean of the observed fluxes (IRodriguez-Pascual et al.l 119971 : lEdelson et al.ll2002l ). 



Finally, column (8) is the ratio of the maximum to minimum flux in the light curves. 



2.4. Time Series Analysis 



A times series analysis of the continuum and H/? light curves was performed to determine 
the mean light travel time lag between continuum and emission line variations. We used two 
cross correlation schemes designed for data sets with uneven time sampling: 



1. An in terpolation scheme (iGaskell &: Sparkelll986l : lGaskell &: Petersonlll987l : IWhite fc Peterson 



19941 ) with an interval of 0.2 day. A cross-correlation function (CCF) is constructed 
from the mean value of the correlation coefficient, r, computed from cross correlating 
both the interpolated line light curve with the original continuum light curve and then 
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the interpolated continuum light curve with the original line light curve, a process 
during which a range of possible lags, r, are imposed on the H/3 light curve. 



2. A time binning scheme (lEdelson fc Kroliklll988l : IWhite &: Petersonlll994l ) with a bin size 
of 1.0 day. Here, a discrete correlation function (DCF) is produced, which determines 
r as a function of lag, similar to the CCF. In this scheme, however, only the actual 
data are cross correlated, and the resulting values of r for all discretely correlated pairs 
are binned as a function of lag. The DCF that results gives the mean value of r in 
each bin, where the corr esponding uncertainty is assigned in a statistical manner (see 



White fc Petersonlll994l ). This method prevents possible spurious lag determinations 



that could potentially arise in the interpolation method due to under-sampling or large 
gaps in the data. 



The resulting CCF and DCF are shown in Figure IH along with the auto-correlation 
function (ACF) computed by cross correlating the continuum with itself. We characterize 
the time delay between the continuum and emission line variations using two parameters 
derived from the CCF; Xpeak is the lag that corresponds to the largest correlation coefficient, 
'"max; ^ud Tcent IS the ceutroid of the CCF based on all points with r > O.Srmax- Time dilation 
corrected values of Tpcak and Tcent determined from the CCF in Figure S] are given in Table 
[3l Uncertainties in both lag parameters are com puted via mode l -indep endent Monte-Carlo 
simulations that employ the bootstrap method of IPeterson et al.l (Il998r ). with the additional 
modifications of IPeterson et al.l ( 200^ 



3. Comparison with Previous Results 

Our measured H/3 lag of Tcent = ^-^^-o'to days from this work is consistent, within the 
errors, to the most recent results for this object by 803, who measured a lag of Tcent = 3.1±1.6 
days. It is not clear how meaningful a direct comparison might be, however, because 803 
measured the time delay between variations in the ~6800A continuum fiux density and the 
integrated Ha fiux rather than between the ~5100A continuum and H/3. We also note that 
the median sampling rate of 803 was larger than our measured lag, suggesting to us that 
the 803 light curves are undersampled. Furthermore, 803 only perform a cross correlation 
analysis based on the DCF method, which sacrifices time resolution. 

Our new time delay measureme nts are inconsist e nt, h owever, with the previous mea- 



surement of Tcent = 5.8l^;8 days by IPeterson et al.l (12004 ) using data from POO. These 



differences are unlikely a luminosity effect, since the average luminosity states of NGC 4051 
were similar during this and the POO campaigns (logALsioo = 41.82 and logALsioo = 41.87, 
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respectively)0. Therefore, we carefully re-examined the light curves used by POO to better 
understand possible causes for the observed inconsistency. 



Netzer fc Maoa (Il990l ) suggested that the cause for a similar inconsi stency between lag 



measurements from two reverberation mapping campaigns of NGC 5548 (INetzer et al.lll990 



Peterson et al.lll99lh was due to different continuum variability timescales observed in the 
separate campaigns: longer continuum variability timescales lead to larger lag measurements. 
However, this explanation is unlikely to be the cause for the current inconsistency between 
our measured lag and that of POO because the prominent variability timescales observed 
in both the POO and current continuum light curves are similar (~40-50 days). Instead, 
the simplest explanation for the inconsistency between our measured lag and that of POO 
is random error. We investigated this possibility by performing Monte Carlo simulations 
using the "Subset 1" H/3 and 5100A continuum flux light curves from POO with the goal 
of estimating the likelihood that a lag of Tcent = 5.8 days would be measured, even if the 



actual BLR radius of the H/5 emission was 2.7 light days, as expected from the iBentz et al. 



( l2009al ) R L relation for the average luminosity of NGC 4051 during this time period. In 
each simulation we created a simulated H/5 light curve by convolving a modified continuum 
light curve with a transfer function that assumed a BLR with a thin spherical shell geometry 
of radius 2.7 light days. The sampling was increased in the modified continuum light curve 
over that of the original Subset 1 continuum light curve by interpolating between the actual 
points on a 0.5 day scale. Noise was added to the flux of each interpolated point using a 
random Gaussian deviate. The size of each deviate was based on the average uncertainty 
in flux of the closest 'real' continuum point on each side of the interpolated point. The 
simulated emission-line light curve was then sampled identically to the original Subset 1 H/3 
light curve. We cross correlated this new emission-line light curve with the original Subset 1 
continuum light curve to determine a reverberation lag. The simulation was repeated 10,000 
times, and lags were measured similarly for each iteration. We found that the average lag 
recovered was 2.7 days (reassuring, since this was our input radius). However, we were 
unable to reproduce even a single lag of 5.8 days. In fact, the largest lag our simulations 
recovered was 4.0 days. A couple of possibihtes suggest themselves: 

• The BLR has physically changed in the 11-year interval between the time of POO's 
Subset 1 and the time of our recent campaign. This is a physical possibility since the 
dynamical timescale of the BLR in NGC 4051 is ^ 5 years. 

• The POO data are undersampled and there are really unresolved variations occurring 



^The average observed flux of tlie 803 campaign was witliin ~ 10% tliat of the Peterson et al. campaign 
as well. 
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on timescales shorter than the typical samphng interval of 2.2 days in Subset 1, the 
best-sampled part of the POO light curve. 

We have no particular reason to believe the former possibility. However, the latter 
is suggested by how POO established the relative uncertainties of their fluxes, namely by 
assuming that there are no true variations on time scales shorter that the typical sampling 
time scales and that any differences between closely spaced measurements reflect random 
errors only, not true variability. The estimates of the relative flux errors in the POO Subset 
1 based on comparing measurements separated by 2 days or less are about 3.2% for both 
the contiuum and the line. In our new data set, obtained with the same instrument, we find 
relative errors of about 1.4% and 2.1% in the continuum and the line, respectively, using the 
method described in Section 12. 3[ We conclude that the flux uncertainties of the POO data 
were overestimated due to short timescale variability. 

Proceeding with the assumption that the POO light curves are undersampled, we isolated 
the portion of the light curves that has the highest sampling across the sharpest features. 
We made this selection in an attempt to avoid occurrences of undersampling more complex 
variability. From the initial light curve, reproduced in Figure \5\ we removed the first 10 
observations that exhibit a broad inflection in the flux with a poorly defined peak. We 
perform a cross correlation analysis on these shortened light curves and determine a shorter 
lag, Tcent = '^■^-I'l days. This lag determination is consistent with both the expected radius 
of 2.7 light days from the R-L relation, the current results, and the results of S03. 

If we continue with the assumption that the light curves of POO are undersampled, 
and the POO flux uncertainties are overestimated, their assigned uncertainties would act to 
decrease the significance of short timescale variability, likely attributing it to noise instead. If 
we impose the average continuum flux uncertainty measured from our current data set (given 
above) on the shortened POO continuum light curve, we can further improve the precision of 
this revised lag measurement to Tcent = '^■^-i^ days. 

We then conducted another simulation in which we applied the sampling rate from the 
POO light curves to the light curves from this work. By undersampling our current light 
curves, we can determine the probability that undersampling could lead to an overestimated 
lag similar to that measured by POO. This type of simulation can provide further evidence 
that the lag measured by POO was an overestimate and a consequence of undersampling. At 
the same time, it could diminish the possibility that the discrepancy in lag measurements 
is due to a difference in the physical conditions or structure of the BLR during the POO 
campaign compared to the present. Using the continuum and H/? light curves shown in Figure 
[3] as the starting point, we modified them similarly to the continuum light curve described for 
our first set of simulations (i.e., increasing the sampling by interpolating between data points 



-14- 



and adding noise to these points with a Gaussian deviate), but this time we interpolated both 
the continuum and emission hne hght curves from this work on a 0.1 day interval We then 
drew sample light curves from this parent light curve with the same length and sampling 
pattern as the full POO light curve shown in Figure O We applied the same cross-correlation 
analysis (as described in Section 12.41) to measure lags from these sample light curves. The 
parent light curves cover a longer time span than the sample light curves and therefore allow 
for multiple iterations of sample light curves to be chosen from different subsets of the parent 
light curves. The first iteration of sample light curves are created from the subset of the 
parent light curves where the beginning points match up, but the ends of the parent light 
curves are discarded. We then build up multiple iterations by shifting the starting point of 
the sample light curve in time by one time step, i.e., 0.1 day, across the parent light curves. 
In this way, we were able to build up 330 sample light curves, where in the last iteration, 
the sample light curves begin in the middle of the parent light curves, but both sets of light 
curves end at the same time. Based on the cross correlation analysis from these 330 sample 
light curves, the probability of measuring Tcent >5.0 days is 0.6% (2 out of 330), and the 
probability of measuring Tcent >3.5 days (i.e, the lag we calculated above from only a portion 
of the POO light curve) is ~8% (25 out of 330). We conclude that undersampling is at least 
a plausible explanation for the difference between the POO results and those reported here. 



4. Black Hole Mass 

Applying virial assumptions to the reverberating gas in the BLR, the mass of the black 
hole can be defined by 



Ml 



BH 



G 



(2) 



where r is the measured emission-line t ime delay, so that c r represents the BLR radius, 
and AV is the BLR velocity dispersion ( jPeterson et al.ll2004l ). The dimensionless factor / 
depends on the structure, kinematics, and inclination of the BLR and is of order unity. 

We estimate the BLR velocity dispersion from the line width of H/3 emission line. The 
line width can be characterized by either the FWHM or the line dispersion, i.e., the sec- 
ond moment of the line profile. The FWHM and the line dispersion, aum, were measured 
from both the mean and the rms spectra of NGC 4051 shown in Figure [H Here, we have 
measured both quant ities and their uncertainties employing methods described in detail by 
Peterson et al.l (120041 ). All measured values of the B.j3 hne width are listed in Tabled Typ- 
ically, the narrow-line emission component of the line should be removed before measuring 
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the line width (see iDenney et aLll2009l ): however, this component could not be reliably iso- 
lated from the rest of the line profile in this object. As a result the line widths measured 
in the mean spectrum, particularly the FWHM, are less reliable for these purposes than the 
widths measured from the rms spectrunj^. 

We calculate the black hole mass for NGC 4051 using Tcent, for the time delay, r, and 
the line dispersion, anne, measured from the H/3 emission line in the rms spectrum, for 
the emission-line width, AV. We utilize the calibration of the reverberation mass scale of 



Onken et al.l (120041 ) for this choice of lag and line width parameters and therefore adopt 
a scale factor value of /=5.5. We then use equation 2 to estimate the black hole mass of 
NGC 4051 to be Mbh = (l-73lo'52) ^ lO^M©. Here, statistical and observational uncertainties 
have been included, but intrinsic uncertainties from sources such as unknown BLR inclination 
cannot be accurately ascertained. 



Marconi et al.l (120081 ) have considered the effect of radiation pressure on SMBH mass 



estimates and provide a new prescription for calculating the SMBH mass that includes a 
correction factor to account for radiation pressure. Radiation pressure acts to partially coun- 
teract the force of gravity on the BLR gas, since the outward radiation force has the same 



radial dependence 



as the inward gravitational force. As a result the BLR gas motions 



are effectively under the influence of a n apparently smaller SMBH mass, which leads to an 
underestimate of the true SMBH mass. iMarconi et al.l (120081 ) determine that although taking 
radiation pressure into account is most important for black holes radiating near the Edding- 
ton limit, it should be considered even for systems with L < L-^dd- On the other hand, in a 
comparis on study of T ype 1 versus Type 2 black hole masses determined with independent 
methods. iNetzerl (120091 ) sees better agreement between the mass and L/L^dd distributions of 
these two populations if radiation pressure forces are neglected. Netzer concludes that either 
the effects of radiation pressure on the BLR gas in these objects is negligible or that BLR 
column densities must be significantly larger, i.e., iVn ^ l O^'^cm"^, than a ssum ed by Marconi 



et al. (2008; = 1 cm ). In a more recent paper, iMarconi et al.l (120091 ) reinvestigate 



the resuhs of iNetzerl j^OoJ and support their findings on the dependence of the effect of 
radiation pressure on column density but conclude that, until it is possible to determine the 
nature of the apparent dependence of A^h on source properties (e.g., L/L^dd), one should 
always consider the possibility that radiation forces are i mportant, and black hole masses 
should consequently be determined using the correction of IMarconi et al.l (120081 ). 



^ Since only BLR emission varies in response to the ionizing continuum on reverberation timescales, flux 
contributions from the narrow-line component will not contaminate the line width measurement in the rms 
spectrum. 
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The importance of radiation pressure forces on black hole mass determinations is still 
under debate. Therefore, in addition to our virial mass estimate for NGC 4051 given above, 
we also estimate t he SMBH mass in N GC 4051 taking radiation pressure into consideration 
(cf. equation 6 of iMarconi et al.ll2008l ). with / = 3.1 ± 1.4 and logg = 7.6 ± 0.3, which are 
derived by Marconi et al. from fits to SMBH masses from reverberation mapping studies. 
With these scale factors and the same line width and BLR radius measurements used above, 
we calculate a mass of MBH-rad = (l-24!l^o;5g) x IO^Mq. Contrary to the expectations from 
the physical arguments, we calculate a mass smaller than that determined in the case where 
we did not consider the effect of radiation pressure. Because NGC 4051 is a low-luminosity 
AGN radiating well below the Eddington limit (L/LEdd = 0.030), the correction to the mass 
due to radiation pressure is small (smaller even than the formal observational uncertainties 
on the mass), adding only 0.26 x IO^Mq to the mass. However, this radiation-corrected mass 
estimate is smaller than our original estimate because th e value of / derived by Marconi et 
al. is a factor of 1.8 smaller than the lOnken et al.l (120041 ) value we adopted above. Because 
this scale factor was derived in a statistical sense by both Onken et al. and Marconi et al., 
the difference in / values between these two studies has the potential to affect the mass of 
an individual object more than would be expected for a statistical sample, particularly for 
the low accretion-rate objects that need only a small radiation pressure correction. 



5. Velocity- Resolved Investigation 



The lag measurements between the continuum and H/3 emission in Table [3] represent the 
average time delay across the BLR. Because the BLR is an extended region and the velocity 
of the gas is most likely a function of position, gas in different locations of the BLR should 
respond to variations in the ionizing continuum flux on slightly different time scales. The 
observable result should be a difference in the reverberation lag measurement in different 
parts of the line profile (i.e., separated in velocity space). Velocity-resolved reverberation 
mapping thus gives us information about the kinematics of gas in the BLR. Previous studies 
of time delay differences between multiple emission lines and the velocity dependence of the 
lag within a single emission line have s hown that the BLR is virialized and commonly contains 



an additional inflow component (e.g.. lGaskelllll988 



1995; Done fc Krohk 1996; Welsh et al 



2007 



Koratkar fc Gaskell 1991; Korista et al. 



Bentz et al.ll2008l ): however, the creation of 



full velocity-delay maps (see iHorne et al.l l2004j ) is an as pect of the reverberation mapping 
technique that has not yet been fully realized (tho ugh see iHorne et al.lll99ll . iDone fc Krolik 
19961 . lUlrich fc Hornd Il996l . and iKoUatschnyl |2003| for previous attempts) . By resolving the 
velocity-dependent reverberation response of the BLR better than has been done in the 
past, we can reconstruct and analyze the velocity-delay map to gain further insights into 
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the geometry and kinematics of the BLR. 

We searched for a velocity-dependent reverberation signal by dividing the H/5 emission 
line flux into 8 velocity-space bins. The blue and red sides of the line were separately 
divided into 4 bins of equal velocity width, covering only the most variable portions of the 
line profile, i.e., the outer-most wavelength boundaries were reduced from those considered 
for the analysis of the full profile to only include flux within the range 4840-4900A (roughly 
±2,000 km s~^). Light curves were created from measurements of the integrated H/3 flux 
in each bin and then cross correlated with the continuum light curve following the same 
procedures described above. The top panel of Figure [6] shows the division of the H/3 line 
profile from the rms spectrum into the eight velocity bins, and the bottom panel shows the 
lag measurements for each of these bins. Error bars in the velocity direction represent the bin 
width. The evidence for a velocity-stratified BLR response to continuum variations is present, 
but marginal. In particular, the lags measured for bins 1-2 and 7-8 are consistent with zero 
and might simply reflect correlated errors due to continuum contamination. Although the 
shape of the velocity-resolved signal in Figure [6] supports our virial assumptions, since the 
higher velocity gas varies on shorter time scales than the low velocity gas, there are no 
strong indications for either outflow or inflow. Outflow could be suggested by the larger 
lag measured in bin 5 (red side of the line) compared to bin 4 (blue side of the line), but 
the difference is very marginal. Even with the high sampling rate and spectral resolution 
we achieved during this campaign, o bserving a v e locity -resolved signal for NGC 4051 as 



clearly as that detected for Arp 151 by lBentz et al.l (120081 ) would have been rather surprising 



for the following reasons. First, the precision with which we can measure a lag is somewhat 
dependent on the median observational sampling interval, which, for NGC 4051, was still not 
much shorter than the measured lag. This indicates that in order to better resolve a velocity- 
dependent signal, we need even higher time resolution for this object. Second, because the 
H/3 line is particularly narrow in this object, there are only a few velocity resolution elements 
across the line. 



6. Discussion 

Based on our simulations, reanalysis of the light curves from POO, and the additional 
arguments we presented above, we conclude that the inconsistency between the past mea- 
surements of the H/3 reverberation lag in NGC 4051 and our present measurements most 
likely results from an overestimation of the lag by POO. Therefore, we adopt the results 
from the current reverberation campaign over previous campaign results measuring this lag 
in NGC 4051. Using the lag measurement of Tcent = ^■87~^o^o days presented here and the 
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TuUy-Fisher distance of iRusselll ( 120031 ) . NGC 4051 is no longer an outlier on the -Rrt,r-I/ 



relation ship. Figure [7] replicates the most recent version of this relationship by lBentz et al. 



(l2009al ) with both the previous and current lag values of NGC 4051 marked. Secure place- 
ment of low-luminosity objects such as NGC 4051 on the i?BLR^-^ relationship is important 
for supporting the extrapolation of this relationship to the even lower-luminosity regime 
potentially populated by intermediate-mass black holes. Additional results from the present 
campaign (Denney et al., in preparat ion), as well as re sults from a recent monitoring cam- 
paign at the Lick Observatory (e.g., iBentz et al.l l2008l ). aim to further populate this low- 
luminosity end of the Rblr^L relationship, thus solidifying the calibration in a relatively 
under-sampled region of the relation. A reliable calibration of this relationship is imperative 
for large studies of black hole masses and galaxy evolution, since it allows for the calculation 
of black hole masses from single-epoch spectra and provides luminosity and radius estimates 
that help constrain parameter space in the search for intermediate-mass black holes. 

Our new results provide a measure of the BLR radius of NGC 4051 that is closer to the 
value naively expected from its luminosity. However, these new results do not resolve the 
unexpected location of this object on the M^yi-L relation: NLSls tend to lie o n a locus of this 



relation with relatively high Eddington ratios (L/LEdd ^ 0.1; see Figure 16 of iPeterson et al. 



20041 ). However, based on the black hole mass of Mbh = (1.731^:^^) X IO^Mq that we have 



calculated, the Eddington ratio of NGC 4051 is still only L/L-^dd = 0.030. It seems unlikely 
that we have overestimated the mass of the black hole (and thus underestimated L/L-^dd), 
since our ma ss measurement is already a factor of a few low er than predicted by the Mbh^ct* 
relationship ( Nelson &: Whittle 1995 : Ferrarese et al.|[200ll ). In this case, the narrowness of 
the Balmer lines in the spectrum of NGC 4051 might be due at least in part to the inclination 
of the BLR relative to our line of sight — indeed, inclination has been invoked as one possible 
way to explain the NLSl phenomenon since t he early days of research on these sources (e.g., 
Osterbrock fc Poggelll985l: I Boiler et al.lll996r) . From observations of the narrow-line region 
in NGC 4051, IChristopoulou et al.l (119971 ) estimate that the inclination of this source is 
~ 50 deg, near the maximum expected for a Type 1 active nucleus in unified models. It is 
entirely reasonable to suppose that the BLR and accretion disk are at the same inclination as 
the narrow-line region. Even so, accounting for this high inclination would increase the line 
width by only about a factor of 2, still within the NLSl regime. Clearly further investigation 
is required to understand the low value of L/L^dd in NGC 4051 compared to other NLSls. 
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Table 1. y-band, Continuum, and H/3 Fluxes for NGC 4051 



JD Fa (5100 A)^ H/5 A4861 Observatory 

(-2450000) (10-15 erg s'^ cm'^ A'^) {lO'^^ erg s'^ cm-^) 



4180.30 


5.03 ± 0.18 




CrAO 


4181.37 


5.66 ± 0.16 




CrAO 


4182.02 


5.47 ± 0.09 




MAGNUM 


4182.41 


5.07 ± 0.22 




CrAO 


4184.79 


5.41 ± 0.20 


5.22 ± 0.11 


MDM 


4185.71 


4.51 ± 0.19 


4.96 ± 0.10 


MDM 


4186.49 


4.51 ± 0.22 




CrAO 


4186.71 


4.89 ± 0.20 


4.15 ± 0.09 


MDM 


4187.38 


4.82 ± 0.22 




CrAO 


4187.85 


5.27 ± 0.20 


4.95 ± 0.10 


MDM 


4188.37 


4.32 ± 0.25 




CrAO 


4188.71 


4.81 ± 0.20 


4.90 ± 0.10 


MDM 


4189.39 


4.75 ± 0.26 




CrAO 


4189.61 


4.78 ± 0.19 


4.62 ± 0.10 


MDM 


4189.96 


5.53 ± 0.21 


4.88 ± 0.10 


MDM 


4189.96 


4.94 ± 0.05 




MAGNUM 


4190.41 


4.82 ± 0.23 




CrAO 


4190.72 


4.82 ± 0.20 


4.73 ± 0.10 


MDM 


4191.38 


4.59 ± 0.32 




CrAO 


4191.62 


4.91 ± 0.20 


4.51 ± 0.09 


MDM 


4191.91 


5.27 ± 0.20 


4.77 ± 0.10 


MDM 


4192.43 


4.84 ± 0.32 




CrAO 


4192.75 


4.02 ± 0.19 


4.44 ± 0.09 


MDM 


4193.61 


4.86 ± 0.20 


4.34 ± 0.09 


MDM 


4193.92 


4.41 ± 0.19 


4.57 ± 0.10 


MDM 


4194.73 


4.62 ± 0.19 


4.60 ± 0.10 


MDM 


4195.45 


4.49 ± 0.16 




UNebr. 


4195.63 


4.30 ± 0.19 


4.06 ± 0.09 


MDM 


4196.62 


4.59 ± 0.19 


4.22 ± 0.09 


MDM 


4197.78 


4.22 ± 0.19 


4.43 ± 0.09 


MDM 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4198.44 


4.34 ± 0.14 




UNebr. 


4198.78 


4.51 ± 0.19 


4.80 ± 0.10 


MDM 


4198.93 


4.46 ± 0.06 




MAGNUM 


4199.36 


4.68 ± 0.17 




CrAO 


4199.40 


4.78 ± 0.15 




UNebr. 


4199.86 


4.55 ± 0.19 


4.37 ± 0.09 


MDM 


4200.40 


4.52 ± 0.15 




CrAO 


4200.72 


4.19 ± 0.19 


4.29 ± 0.09 


MDM 


4201.31 


4.84 ± 0.23 




CrAO 


4201.73 


4.59 ± 0.19 


4.49 ± 0.09 


MDM 


4202.37 


4.56 ± 0.18 




CrAO 


4202.95 


4.93 ± 0.08 




MAGNUM 


4204.40 


4.53 ± 0.17 




CrAO 


4204.73 


4.35 ± 0.19 


4.58 ± 0.10 


MDM 


4205.34 


4.34 ± 0.16 




CrAO 


4205.62 


3.97 ± 0.18 


4.37 ± 0.09 


MDM 


4205.83 


4.32 ± 0.03 




MAGNUM 


4205.91 


4.28 ± 0.19 


4.54 ± 0.10 


MDM 


4206.36 


4.19 ± 0.18 




CrAO 


4206.40 


4.23 ± 0.23 




UNebr. 


4206.62 


4.31 ± 0.19 


4.46 ± 0.09 


MDM 


4207.40 


4.30 ± 0.15 




UNebr. 


4207.82 


4.40 ± 0.19 


4.60 ± 0.10 


MDM 


4208.34 


4.42 ± 0.15 




CrAO 


4208.40 


4.26 ± 0.17 




UNebr. 


4208.62 


4.24 ± 0.19 


4.58 ± 0.10 


MDM 


4208.88 


4.25 ± 0.03 




MAGNUM 


4209.40 


4.46 ± 0.15 




CrAO 


4209.78 


4.81 ± 0.20 


4.81 ± 0.10 


MDM 


4210.62 


5.01 ± 0.20 


4.99 ± 0.10 


MDM 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4211.41 


4.70 ± 0.31 




CrAO 


4212.34 


4.64 ± 0.15 




CrAO 


4212.72 


5.11 ± 0.20 


5.01 ± 0.11 


MDM 


4212.75 


4.78 ± 0.03 




MAGNUM 


4213.34 


5.11 ± 0.19 




CrAO 


4213.74 


4.90 ± 0.20 


4.50 ± 0.09 


MDM 


4214.34 


4.65 ± 0.20 




CrAO 


4214.73 


4.92 ± 0.20 


5.15 ± 0.11 


MDM 


4215.40 


4.72 ± 0.21 




CrAO 


4215.74 


4.84 ± 0.20 


4.84 ± 0.10 


MDM 


4216.32 


5.15 ± 0.21 




CrAO 


4216.73 


4.73 ± 0.19 


5.09 ± 0.11 


MDM 


4217.35 


4.49 ± 0.23 




CrAO 


4217.73 


4.73 ± 0.19 


5.22 ± 0.11 


MDM 


4218.31 


4.58 ± 0.19 




CrAO 


4218.80 


4.31 ± 0.19 


4.74 ± 0.10 


MDM 


4218.91 


4.70 ± 0.04 




MAGNUM 


4219.32 


4.73 ± 0.21 




CrAO 


4219.40 


4.71 ± 0.16 




UNebr. 


4219.83 


5.06 ± 0.20 


4.70 ± 0.10 


MDM 


4220.31 


4.78 ± 0.25 




CrAO 


4220.40 


5.05 ± 0.16 




UNebr. 


4220.74 


5.18 ± 0.20 


4.89 ± 0.10 


MDM 


4221.35 


5.23 ± 0.38 




CrAO 


4221.74 


4.91 ± 0.20 


4.75 ± 0.10 


MDM 


4221.99 


4.78 ± 0.08 




MAGNUM 


4222.40 


5.08 ± 0.26 




CrAO 


4222.74 


4.91 ± 0.20 


4.83 ± 0.10 


MDM 


4223.38 


4.47 ± 0.25 




CrAO 


4223.74 


4.73 ± 0.19 


5.22 ± 0.11 


MDM 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4224.37 


5.12 ± 0.19 




CrAO 


4224.73 


4.81 ± 0.20 


5.41 ± 0.11 


MDM 


4225.35 


4.46 ± 0.16 




CrAO 


4225.80 


4.30 ± 0.19 


4.56 ± 0.10 


MDM 


4225.88 


4.43 ± 0.03 




MAGNUM 


4226.29 


4.46 ± 0.20 




CrAO 


4226.75 


3.88 ± 0.18 


4.55 ± 0.10 


MDM 


4227.43 


4.28 ± 0.23 




CrAO 


4227.74 


4.12 ± 0.19 


4.58 ± 0.10 


MDM 


4228.80 


3.81 ± 0.18 


4.54 ± 0.10 


MDM 


4229.37 


3.91 ± 0.15 




CrAO 


4229.78 


3.91 ± 0.18 


4.31 ± 0.09 


MDM 


4230.73 


4.38 ± 0.19 


4.36 ± 0.09 


MDM 


4231.36 


4.07 d= 0.16 




CrAO 


4231.59 


4.58 ± 0.28 




UNebr. 


4231.75 


4.81 ± 0.20 


4.61 ± 0.10 


MDM 


4232.28 


4.47 ± 0.16 




CrAO 


4232.38 


4.59 ± 0.22 




UNebr. 


4232.73 


4.22 ± 0.19 


4.53 ± 0.09 


MDM 


4233.32 


4.51 ± 0.15 




CrAO 


4233.44 


4.88 ± 0.18 




UNebr. 


4233.73 


4.64 ± 0.19 


4.72 ± 0.10 


MDM 


4234.32 


4.17 ± 0.16 




CrAO 


4234.73 


4.51 ± 0.19 


4.69 ± 0.10 


MDM 


4234.85 


4.29 ± 0.03 




MAGNUM 


4235.31 


4.35 ± 0.14 




CrAO 


4235.46 


4.80 ± 0.36 




UNebr. 


4235.73 


4.43 ± 0.19 


4.56 ± 0.10 


MDM 


4236.31 


4.53 ± 0.15 




CrAO 


4236.73 


3.96 ± 0.18 


4.54 ± 0.09 


MDM 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4237.31 


4.17 ± 0.14 




CrAO 


4237.52 


4.13 ± 0.28 




UNebr. 


4237.73 


3.94 ± 0.18 


4.28 ± 0.09 


MDM 


4238.49 


4.38 d= 0.17 




UNebr. 


4238.73 


4.16 ± 0.19 


4.42 ± 0.09 


MDM 


4238.93 


4.05 ± 0.06 




MAGNUM 


4239.35 


4.19 ± 0.17 




CrAO 


4239.43 


4.34 ± 0.18 




UNebr. 


4239.75 


3.81 ± 0.18 


4.40 ± 0.09 


MDM 


4240.29 


3.69 ± 0.14 




CrAO 


4240.47 


3.69 ± 0.18 




UNebr. 


4240.72 


3.97 ± 0.18 


4.37 ± 0.09 


MDM 


4241.31 


3.74 ± 0.16 




CrAO 


4241.38 


3.92 ± 0.24 




UNebr. 


4241.73 


3.75 ± 0.18 


4.30 ± 0.09 


MDM 


4242.29 


3.75 ± 0.14 




CrAO 


4242.75 


3.56 ± 0.18 


3.91 ± 0.08 


MDM 


4243.29 


3.35 ± 0.19 




CrAO 


4243.74 


3.39 ± 0.18 


4.05 ± 0.09 


MDM 


4244.78 


3.68 ± 0.18 


3.90 ± 0.08 


MDM 


4245.35 


4.38 ± 0.25 




CrAO 


4245.75 


4.36 ± 0.19 


4.03 ± 0.09 


MDM 


4245.77 


4.21 ± 0.06 




MAGNUM 


4246.34 


4.95 ± 0.18 




CrAO 


4246.40 


4.49 ± 0.17 




UNebr. 


4246.74 


4.06 ± 0.19 


4.12 ± 0.09 


MDM 


4247.74 


4.57 ± 0.19 


4.20 ± 0.09 


MDM 


4248.33 


4.50 ± 0.21 




CrAO 


4248.73 


4.50 ± 0.19 


4.48 ± 0.09 


MDM 


4249.45 


4.36 ± 0.22 




CrAO 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4249.74 


4.40 ± 0.19 


4.52 ± 0.09 


MDM 


4250.41 


3.56 ± 0.33 




CrAO 


4250.74 


4.07 ± 0.19 


4.44 ± 0.09 


MDM 


4251.32 


3.74 ± 0.30 




CrAO 


4251.74 


4.02 ± 0.19 


4.39 ± 0.09 


MDM 


4252.39 


4.47 ± 0.27 




CrAO 


4252.73 


4.15 ± 0.19 


4.22 ± 0.09 


MDM 


4252.88 


4.17 ± 0.05 




MAGNUM 


4253.73 


4.44 ± 0.19 


4.19 ± 0.09 


MDM 


4254.35 


4.62 ± 0.16 




CrAO 


4254.39 


4.52 ± 0.17 




UNebr. 


4255.38 


4.73 ± 0.24 




CrAO 


4255.76 


4.40 ± 0.19 


4.50 ± 0.09 


MDM 


4256.35 


5.02 ± 0.14 




CrAO 


4256.71 


4.75 ± 0.19 


4.56 ± 0.10 


MDM 


4257.40 


4.80 ± 0.17 




CrAO 


4257.74 


4.64 ± 0.19 


4.47 ± 0.09 


MDM 


4258.35 


4.52 ± 0.16 




CrAO 


4258.50 


5.06 ± 0.32 




UNebr. 


4258.76 


4.83 ± 0.20 


4.91 ± 0.10 


MDM 


4259.34 


5.15 ± 0.17 




CrAO 


4259.75 


4.92 ± 0.20 


4.97 ± 0.10 


MDM 


4259.84 


4.81 ± 0.13 




MAGNUM 


4260.30 


4.96 ± 0.17 




CrAO 


4260.75 


4.36 ± 0.19 


4.70 ± 0.10 


MDM 


4261.31 


5.08 ± 0.16 




CrAO 


4261.42 


5.15 ± 0.14 




UNebr. 


4261.74 


4.82 ± 0.20 


5.16 ± 0.11 


MDM 


4262.30 


4.78 ± 0.16 




CrAO 


4262.74 


4.39 ± 0.19 


4.72 ± 0.10 


MDM 
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Table 1 — Continued 



JD 
-2450000) 


Fa (5100 A)^ 
(10"^^ erg cm~^ A~^) 


H/5 A4861 
(10~^^ erg cm~^) 


Observatory 


4263.35 


4.68 ± 0.17 




CrAO 


4263.39 


4.47 ± 0.22 




UNebr. 


4263.72 


4.52 ± 0.19 


5.01 ± 0.11 


MDM 


4263.88 


4.61 ± 0.11 




MAGNUM 


4264.76 


4.73 ± 0.19 


5.12 ± 0.11 


MDM 


4265.76 


5.43 ± 0.20 


5.15 ± 0.11 


MDM 


4266.34 


4.84 ± 0.27 


5.18 ± 0.23 


CrAO 


4266.76 


5.55 ± 0.21 


5.29 ± 0.11 


MDM 


4267.30 


4.43 ± 0.26 


5.43 ± 0.24 


CrAO 


4267.75 


4.55 ± 0.19 


5.24 ± 0.11 


MDM 


4268.29 


4.33 ± 0.26 


5.34 ± 0.23 


CrAO 


4268.75 


4.38 ± 0.19 


5.17 ± 0.11 


MDM 


4269.32 


4.57 ± 0.26 


5.16 ± 0.23 


CrAO 


4269.75 


4.56 ± 0.19 


5.39 ± 0.11 


MDM 


4269.84 


4.37 ± 0.06 




MAGNUM 


4270.36 


5.09 ± 0.27 


5.12 ± 0.23 


CrAO 


4270.37 


4.37 ± 0.20 




UNebr. 


4271.31 


4.52 ± 0.26 


5.38 ± 0.24 


CrAO 


4272.85 


4.54 ± 0.04 




MAGNUM 


4274.31 


4.78 ± 0.26 


5.16 ± 0.23 


CrAO 


4275.81 


4.31 ± 0.03 




MAGNUM 


4276.41 


4.80 ± 0.20 




UNebr . 


4277.29 


4.42 ± 0.26 


4.79 ± 0.21 


CrAO 


4278.29 


4.01 ± 0.25 


4.99 ± 0.22 


CrAO 


4278.48 


3.66 ± 0.37 




UNebr. 


4278.84 


4.14 ± 0.06 




MAGNUM 


4279.31 


4.33 ± 0.26 


4.74 ± 0.21 


CrAO 


4279.32 


4.22 ± 0.21 




CrAO 


4280.32 


4.61 ± 0.26 


4.87 ± 0.21 


CrAO 


4280.34 


4.38 ± 0.43 




CrAO 
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Table 1 — Continued 



JD 

(— 245UUUUj 


b\ (5100 A) 

/'in— 15 „ —1 —2 A — 1\ 

(10 erg s cm A j 


TT o \ /I or* 1 

Hp A4861 

/'in— 13 „ —1 —2\ 

(10 erg s cm ) 


Observatory 


4281 32 


4 34 ± 24 




CrAO 


4281.32 


4 38 ± 26 


4.88 ± 0.22 


CrAO 


4282.33 


3.96 ± 0.24 




CrAO 


4282.35 


4.22 ± 0.25 


4.47 ± 0.20 


CrAO 


4283.31 


3 78 ± 25 


4 63 ± 20 


CrAO 


4283 33 


3 75 ± 24 

• 1 — 1 — vy • -L. 




CrAO 


4283.42 


3 56 ± 20 




UNebr. 


4284.30 


3.64 ± 0.24 


4 50 ± 20 


CrAO 


4284.31 


3 87 ± 18 




CrAO 


4285.77 


3 70 ± 09 




MAGNUM 


4289.29 


4.20 ± 0.25 


4.22 ± 0.19 


CrAO 


4289.45 


3 65 ± 25 




UNebr. 


4290 30 


3 85 ± 25 


4.41 ± 0.19 


CrAO 


4290.40 


4.27 ± 0.22 




UNebr. 


4291 30 


3 53 ± 24 

— 1 — vy«^^ 


4 34 ± 19 


CrAO 


4294.29 


4.63 ± 0.18 




CrAO 


4295.32 


4.51 ± 0.34 




CrAO 


4296.29 


4.35 ± 0.26 


4.67 ± 0.21 


CrAO 


4296.30 


4.57 ± 0.17 




CrAO 


4298.32 


4.78 ± 0.26 


4.38 ± 0.19 


CrAO 


4299.28 


4.29 ± 0.26 


4.65 ± 0.20 


CrAO 


4299.31 


4.87 ± 0.16 




CrAO 


4300.28 


4.14 ± 0.25 


4.68 ± 0.21 


CrAO 


4304.79 


4.66 ± 0.11 




MAGNUM 


4311.76 


4.79 ± 0.09 




MAGNUM 



'^This column contains the average continuum flux density measured at rest- 
frame ~5100A from spectroscopic observations as well as the F-band flux from 
photometric observations. Spectroscopic and photometric fluxes were intercali- 
brated and merged to create a single continuum light curve (see Section 2.3). 



-33- 







Table 2. 


Light Curve Statistics 










Sampling 




Mean 






Time 




Interval(days) 


Mean 


Fractional 






Series 


N 


(T^ T'median 


Flux^ 


Error 


^ vax 




(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


5100 A 


186 


0.71 0.56 


4.5 ±0.4 


0.04 


0.09 


1.69 ±0.11 


R/3 


100 


1.17 1.00 


4.7 ±0.3 


0.02 


0.07 


1.39 ±0.04 



Same flux units as Table 1 for 5100 A continuum and H/3, respectively. 
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Table 3. Reverberation Results 



Parameter 


Value 


(1) 


(2) 



'''cent 
''"peak 

criine(mean) 
FWHM (mean) 

criine(rms) 
FWHM frms) 



'BH 



Ml 



BH-rad 



1.87+°J^ days 



2.60+^;^;; days 



0.79 
1.40 

1045 ± 4 km s'^ 
799 ± 2 km s^^ 
927 ± 64 km s^^ 
1034 ±41 km s'^ 
:i.73l°:|) X lO^Me 
:i.24+°-^I) X lO^M^ 



'^Using lOnken et al.l (120041 ) calibra- 
tion. 



Using iMarconi et al.l ( 20081 ) calibra- 
tion. 
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Fig. 1. — Mean and rms spectra of NGC 4051 from MDM observations. The rms spectrum 
was formed after removing the [O ill] A4959 and [O ill] A5007 narrow emission hnes. The 
variabihty signature of H/5 is clearly visible in the rms spectrum, and the large increase in 
rms flux shortward of 4800 A is due to variations in the broad Hell A4686 emission line. 
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Fig. 2. — Light curves showing complete set of observations from all four sources. The 
top panel shows the 5100 A continuum flux in units of 10~^^ erg s~^ cm~^ while the 

bottom is the H/5 A4861 line flux in units of 10~^^ erg cm~^. The open triangles (CrAOsp) 
correspond to spectroscopic observations taken at CrAO, while the closed triangles (CrAOph) 
represent photometric observations from CrAO. 
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Fig. 3. — Same as Fig. [2] except data from all sources have been merged and closely spaced 
observations binned such that weighted averages were calculated for continuum observations 
separated by less than 0.25 day and H/3 observations separated by less than 0.5 day. 
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-20 -10 10 20 



T(days) 

Fig. 4. — Cross correlation function (CCF; solid line), discrete correlation function (DCF; 
filled circles), and auto correlation function (ACF; dotted line) from time series analysis of 
the continuum and H/3 light curves of NGC 4051 shown in Fig. [3l 
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Fig. 5. — Optical continuum and H/3 light curves reproduced from POO. Points with X's 
represent those that were excluded from the new time series analysis of this light curve 
described in Section [31 Units are the same as in Fig. [2l 
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Fig. 6. — Velocity-resolved B.(3 rms spectrum profile (top) and time-delay measurements 
(bottom) for NGC 4051. Vertical dashed lines plotted on the line profile (top) and error bars 
on the lag measurements in the velocity direction (bottom) show the bin size, with each bin 
labeled by number in the top panel. Error bars on the lag measurements are determined 
similarly to those for the mean BLR lag. The horizontal solid and dotted lines in the bottom 
panel show the mean BLR centroid lag and associated errors, calculated in Section [231 while 
the horizontal dotted-dashed line in the top panel represents the linearly-fit continuum level. 
Flux units are the same as in Fig. [TJ 
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Fig. 7. — Most recently calibrated Rblr^L relation f Bentz et al.1 2009a. solid line). The filled 
square shows the location of NGC 4051 based on results from lPeterson et al.l (120041 ) and used 
by Bentz et al. The filled circle shows the new lag measurement of 1.87 days presente d in this 
work at the luminosity calculated using the Tully-Fisher distance to NGC 4051 of iRussell 
( 120031 ): the error bar in luminosity reflects both the range of flux variations, as for all of the 
other data points, plus the uncertainty due to the distance, added in quadrature. The error 
bar i s asymmetri c , as w' e favor the larger distance. Open squares represent other objects 
from iBentz et al.1 (l2009ah . 



